Spatial Pauli-blocking of spontaneous emission in optical lattices 
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Spontaneous emission by an excited fermionic atom can be suppressed due to the Pauli exclusion 
principle if the relevant final states after the decay are already occupied by identical atoms in the 
ground state. Here we discuss a setup where a single atom is prepared in the first excited state 
on a single site of an optical lattice under conditions of very tight trapping. We investigate these 
phenomena in the context of two experimental realizations: (1) with alkali atoms, where the decay 
rate of the excited state is large and (2) with alkaline earth-like atoms, where the decay rate from 
metastable states can be tuned in experiments. This phenomenon has potential applications towards 
reservoir engineering and dissipative many-body state preparation in an optical lattice. 



PACS numbers: 03.75.Ss, 37.10.Jk, 42.50. Ct, 67.85.-d 



I. INTRODUCTION 

An atom in free space prepared in an excited electronic 
state will decay by spontaneous emission to lower lying 
electronic states. Spontaneous emission is a basic ingredi- 
ent in light scattering from atoms, and processes involv- 
ing spontaneous emission constitute fundamental decay, 
and thus decoherence mechanisms, when manipulating 
atoms with laser light. There is an extensive literature 
on suppression of the (free space) spontaneous emission 
rate. This can be achieved first of all by engineering 
the density of states of the radiation modes of the emit- 
ted photons, e.g. by placing emitters in a cavity [T], in 
a photonic bandgap material |^, or close to a surface 
[3]. On the other hand, recent experimental advances 
with quantum degenerate Fermi gases have motivated 
theoretical studies of "Pauli-blocking" spontaneous emis- 
sion, where the spontaneous emission of an atom from 
the first excited electronic state to the ground state is 
blocked by other fermionic atoms occupying possible fi- 
nal states into which the atom can decay. This scenario 
and prospects for experimental observation have so far 
been discussed for a single excited atom above a Fermi 
sea of many trapped ground state atoms [H |5], where 
the condition for Pauli-blocking spontaneous emission is 
Ep > En with Ep the Fermi energy of ground state 
atoms, and Er — h^k'^/2M the recoil energy, where k 
denotes the wavenumber of the emitted photons and M 
the atomic mass. 

Here we will instead discuss a complementary and con- 
ceptually simpler setup, in which a single atom is pre- 
pared in the first excited state on a single site of an optical 
lattice under conditions of very tight trapping realizing 
the Lamb Dicke limit, i.e. ry = 2t:x()/Xl ^ 1 with xq 
the size of the atomic wavepacket and A^ the wavelength 
of the trapping laser. In this case under suitable con- 
ditions, which we discuss below, even a single fermionic 
atom prepared on the same lattice site in the electronic 
ground state and motional ground state matching the 



atom in the excited state can block the dominant decay 
channel for the excited state to the ground state. Below 
we will analyze this suppression of decay, in particular 
in light of the new opportunities opened by recent ex- 
periments and theoretical proposals involving fermionic 
alkaline earth atoms [5HT5] , but also with alkali atoms in 
optical lattices [TiHTO] . 

The paper is organized as follows. In Sec. |ll] we give 
a brief qualitative overview of the atomic decay dynam- 
ics and introduce two experimental schemes to observe 
Pauli-blocked spontaneous emission. Details of these sce- 
narios are discussed in Sec. |III[ In Sec. |IV] we study the 
properties of the emitted photon wavepacket and the pos- 
sibility to shape its form by preparing the blocking atom 
in appropriate initial states. We conclude with a sum- 
mary in Sec. [viand also give an outlook towards possible 
applications of Pauli-blocked spontaneous emission as a 
tool for reservoir engineering in the context of dissipative 
preparation of many-body quantum states and phases of 
fermions in optical lattices [T7H25] . 



II. OVERVIEW 

In this section, we first give a qualitative description of 
the effect of Pauli-blocked spontaneous emission in opti- 
cal lattices and subsequently outline possible experimen- 
tal realizations using alkali and alkaline-earth atoms. 

Let us consider a three dimensional optical lattice that 
is loaded with two identical fermionic atoms per site and 
which is sufficiently deep that tunneling between neigh- 
boring sites is negligible on experimental timescales. Let 
us further assume that on each site one of the atoms is in 
an electronically excited state |e) while the other resides 
in its internal ground state \g). Due to its coupling to the 
electromagnetic field, the excited atom can undergo ra- 
diative decay. As a result of the Pauli exclusion principle, 
the presence of the ground state atom reduces the num- 
ber of available decay channels, prolonging the lifetime 
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FIG. 1. (a) Feff ^ v: The two atoms prepared in the internal 
excited state |e) and the internal ground state \g) feel the 
same external harmonic oscillator potential. In the Lamb- 
Dicke limit (?7 <C 1), the dominant decay channel for a single 
particle |e) |0) — >■ Ig) |0) with a rate To is now blocked by the 
ground state atom due to the Pauli exclusion principle. The 
excited atom decays under change of its motional state with 
a rate Fefi that is reduced of order rf . (b) Feff S> v. In this 
scenario, the notion of motional eigenstates is meaningless for 
atoms in |e). Yet again, both atoms can be prepared in the 
same initial motional state represented by the wavepacket jO) 
(solid line). Dispersion of the motional wavepacket for the 
excited atom (dashed line) is slow compared to the effective 
decay dynamics. As before, the dominant decay channel is 
blocked and the effective decay rate is reduced of order rf . 



of the excited state. For this blocking effect to be signifi- 
cant we require that both atoms are initially prepared in 
the same motional state, that the atoms are tightly con- 
fined to a region smaller than the optical wavelength of 
the electronic transition (Lamb-Dicke regime) and that 
decay from the excited state |e) is only possible to the 
internal ground state \g) (see Fig. [TJa)). 

In the regime without tunneling the spatial potential 
on each lattice site is well approximated by a three dimen- 
sional harmonic oscillator potential. Here, for simplicity 
of this qualitative discussion, we regard motional exci- 
tations only in one dimension where the oscillator trap- 
ping frequency is v and the vibrational eigenstates are 
denoted \n) {n = 0, 1, . . .). This will later be extended 
to a 3D model. Furthermore, we assume for the moment 
that atoms in the internal excited state feel the same 
oscillator potential as ground state atoms, and that the 
natural linewidth F of the state |e) is small compared to 
u. While these requirements simplify the discussion they 
are not mandatory, as we will argue below. For our initial 
state of interest, which is denoted |?/;o) = CggCeo l^ac) in 
second quantized notation, both atoms reside in the mo- 
tional ground state |0) and one is internally excited (see 
Fig. flla)). Here, cj„ (Ce„) creates a particle on the site 
we consider in the electronic state l^) (|e)) and motional 
state \n) . Treating the system within a Weisskopf-Wigner 
approximation |26| . the initial state |V'o) couples to the 



possible final atomic states Cgoc],„ 



vac) under emission of 



a one-photon wavepacket. The resulting dynamics within 
this ansatz yields an exponential decay for the population 
in the excited state with an effective decay rate 

roff = r^/' duN{u)\{n\e-'-^('^+^^')\0)\\ (1) 

In Eq. ([l]), F = d^gLOQ/{iT:£Qhc^) is the usual single par- 
ticle decay rate from |e) to \g) where deg and wq are 
the dipole matrix element and frequency of the optical 
transition, respectively; a^ (a) denotes the creation (an- 
nihilation) operator for the harmonic oscillator and N{u) 
is the angular distribution of dipole radiation projected 
onto the oscillator axis. The parameter f] = 2'kxq/\q 
is given by the ratio betw een the oscillator ground state 
length a;o = \Jh/{2Mi') with atomic mass M and the 
wavelength of the emitted light Aq — I-kcJujq. Note that 
in the sum over the final vibrational modes for the de- 
caying atom, the mode n = is excluded due to the 
fermionic property cJ.qc],q = 0, which reflects the Pauli- 
blocking of this particular channel. 

In the Lamb-Dicke regime xq <^ Aq the exponential in 
Eq. (IT]) can be expanded in the small parameter 77 <C 1. 
In this limit, the leading contribution results from decay 
into the vibrational mode n = 1, yielding an effective de- 
cay rate TcS ~ r/^F. In contrast, a single particle decays 
at a total rate F. For two atoms, the excited state's de- 
cay rate is decreased by a factor of order 77^ . This results 
from blocking of the dominant decay channel, which cor- 
responds to a decay while preserving the motional state, 
due to Pauli's exclusion principle. 

Let us now mention details and additional effects that 
modify the simplified picture outlined so far. Coming 
back to the preparation of the initial state c^gCgQ |vac), 
we distinguish the two scenarios Feff <C i" and Feff ^ v. 
For Feff <C ly (see Fig. [ija)) the decay is slow compared 
to the atomic motion. Such a situation can be realized by 
using metastable excited states of alkaline-earth atoms, 
where it is also possible to produce equal trapping poten- 
tials for ground state and excited atoms via magic wave- 
length lattices [H [57] . In this sideband resolved regime 
it is possible to laser-excite one of the atoms while ex- 
clusively coupling to the motional ground state. The 
opposite regime Feg » v (see Fig. ^h)) is naturally 
realized for short-lived low-lying excited states in alkali 
atoms. Here, laser excitation from \g) to |e) is performed 
in the strong-excitation regime |28) where the initial mo- 
tional wavepacket remains essentially unchanged during 
the electronic transition and corresponds to the motional 
ground state for atoms in the internal state \g). In this 
regime a possibly different (anti-) trapping potential for 
excited state atoms gives rise only to a small correction 
to the Pauli-blocked decay dynamics. In this way, it is 
possible to initially prepare the ground and excited atom 
in the same motional state for either scenario, which will 
therefore lead to an effective decay rate Fgff ^ •q'^T in 
either case. 

Furthermore, additional two-body effects appear. 
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FIG. 2. Relevant level structure (not to scale) and experi- 
mental sequence for alkaline earth-like atoms (^'^^Yb in this 
example). (I) Two atoms per lattice site are prepared in \g t) 
and \g \), respectively; (II) one atom is excited from \g \.) to 
the metastable state |e f) under a flip of its nuclear spin; (III) 
spontaneous emission from |e t) is induced by admixing the 
fast decaying state |J t)- Nuclear spin flips during the induced 
decay are suppressed by applying a large external magnetic 
fleld which decouples nuclear and electronic spin in the ^Pi 
manifold. Consequently the induced decay from |e t) to \g t) 
is Pauli-blocked. 



which will be taken into account in a more accurate treat- 
ment based on a second quantized many particle master 
equation (Section |III[ ). Since both atoms are confined to 
a region much smaller than Aq, their mutual coupling to 
the radiation field gives rise to dipole-dipole interaction 
between the two atoms. In addition, cross-damping pro- 
cesses, which in other contexts are also responsible for 
super- and subradiance, can modify the decay character- 
istics. These effects are more often encountered in the 
limit F ^ i^. Additionally, the influence of state chang- 
ing collisions between the atoms will be addressed. 

In the following two subsections we give an overview of 
possible realizations of Pauli-blocked spontaneous emis- 
sion with ultracold atoms in optical lattices using (A) 
alkaline earth atoms and (B) alkali atoms. We will dis- 
cuss the experimental requirements and steps for each 
scenario and show how the electronic dynamics of such 
multi-level atoms can be related to the effective two-level 
scheme discussed above. 



A. Alkaline earth atoms 

Alkaline earth atoms with two valence electrons of- 
fer the opportunity to demonstrate the effect of Pauli- 
blocked spontaneous emission in a way that is closely 
related to the qualitative model outlined above. These 
atoms have singlet ground states and metastable triplet 
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FIG. 3. Relevant level structure (not to scale) and experi- 
mental sequence for alkali atoms (^''K in this example). (I) 
Two atoms per lattice site are prepared in \g) and \i), respec- 
tively; (II) flrst, one atom is excited from \g) to the excited 
state |e); (III) then the other atom is transferred from \i} to 
\g). This pulse sequence is fast on the timescale given by F"^. 
Choosing states with maximal mp values ensures that \g) and 
|e) form a closed two-level system. The decay of the excited 
atom back to \g) is Pauli-blocked. 



states with lifetimes on the order of several tens of 
seconds. Since this exceeds accessible experimental 
timescales, we will use such a state as the excited state 
in our scheme and induce spontaneous emission by weak 
laser coupling to an intermediate state which rapidly de- 
cays to the internal ground state. Hence, the effective 
tunable decay rate can be made significantly larger than 
the natural linewidth of the metastable state. By oper- 
ating the optical lattice at the "magic" wavelength it is 
possible to have equal trapping potentials for the ground 
and metastable excited state [Hl[57]. Furthermore, the 
fermionic isotopes have an additional nuclear spin degree 
of freedom, which we will make use of in the context of 
initial state preparation. 

The relevant electronic level structure and experimen- 
tal sequence I through III is illustrated in Fig. [2] (I) Af- 
ter the lattice is adiabatically ramped up on an ultracold 
cloud of atoms, two atoms per lattice site with the nu- 
clear spin projections |t) and \D reside in their motional 
and internal ground state \g) |0) forming a band insula- 
tor. (II) With a TT pulse on the weakly dipole-allowed 
transition \g) — > |e), one of the atoms is transferred from 
\g I) |0) to the excited state |e f) |0) while flipping its nu- 
clear spin. (Ill) Finally, spontaneous emission from |e) is 
induced by admixing a fast decaying intermediate state. 
In order to avoid that the excited atom decays via the 
channel |e f) — > 1.9 i) which is not Pauli-blocked, one can 
either use a sufficiently large external magnetic field [29] , 
or alternatively choose an appropriate polarization of the 
dressing laser. For both cases this leads to a decoupling 
of electronic and nuclear spin for the induced decay pro- 
cess (see Fig. [2]), and as the intermediate state is only 
virtually populated, |e f) and \g f) form a closed system 
regarding the internal degrees of freedom. 

As a consequence the described system constitutes a 
realization of the model system illustrated in Fig. flla) 



where the excited atom can only undergo decay under 
change of its motional state. The induced decay rate is 
by a factor of order rf' smaller than it is for a single atom. 



B. Alkali atoms 

In contrast to the alkaline earth case, alkali atoms have 
one valence electron whose lower excited states decay 
to the ground state on dip ole- allowed transitions with 
a large linewidth F on the order of several MHz. 

In Fig. |3] a typical alkali level structure and a possible 
experimental scheme is shown. We choose to make use of 
the transition between the hyperfine sublevels |e) and \g) 
with maximal mp values, because the only dipole-allowed 
transition for atoms prepared in |e) is decay to \g), so 
these states form an effective closed two-level system. (I) 
The lattice is prepared with two atoms per lattice site 
in the ground state sublevels |i) and \g) (both in their 
motional ground state |0)). (II) The first atom is excited 
with a state selective tt pulse from \g) to |e). (Ill) The 
second atom is transferred from |i) to \g) with a radio 
frequency or Raman tt pulse. The total pulse sequence 
should be fast compared to the timescale F^^ so that 
spontaneous emission during the preparation sequence 
can be neglected. The excited atom then undergoes ra- 
diative decay with a Pauli-blocked rate on the order of 
r?2F. 

As a consequence of the large linewidth F imperfec- 
tions are introduced that will qualitatively modify the 
Pauli-blocking effect: (i) the induced dipole-dipole in- 
teraction between atoms is strong for large F and will 
excite motional states that are not blocked; (ii) due to 
the fast excitation pulses the transfer to |e) is performed 
in the strong-excitation regime where the Rabi frequency 
is much larger than the oscillator spacing v of the ground 
state. Thus, the initial motional state |0) is displaced by 
the momentum recoil of the absorbed photon. This leads 
to slightly different motional states for the ground and ex- 
cited state atom; (iii) in general, the atoms encounter dif- 
ferent (anti-)trapping potentials for their internal ground 
and excited states. However, due to the large separation 
of timescales i^ ^ F, corrections to the dynamics due 
to atomic motion during the decay are small. All these 
effects will be discussed in more detail in the following 
section. 



III. ANALYSIS AND RESULTS 

In this section we first discuss the many-body master 
equation for fcrmionic two-level atoms in second quanti- 
zation. We show how Pauli-blocking emerges naturally in 
this description and discuss collective effects which arise 
for more than one atom. Subsequently, we present a de- 
tailed quantitative discussion of the experimental realiza- 
tions with both alkaline earth-like and alkali atoms. 



A. Master equation 

We consider fcrmionic two-level atoms with internal 
states |e) and \g) separated by the energy /iwo = hck^ 
and coupled to the radiation field acting as a bath. The 
center of mass motion of the atoms takes place in a three 
dimensional harmonic oscillator potential with trapping 
frequencies u = (i^i, i^2j I'a) foi' the three spatial dimen- 
sions and eigenstates \n) = \ni) \n2) \n^) {uj = 0, 1, • • • ) 
on each lattice site for a deep optical lattice with strongly 
suppressed tunneling. We are interested in a situation 
where the atoms are confined to a region small compared 
to the optical wavelength Aq — 27r/fco. Under assumption 
of the Born-Markov approximation, the dynamics is de- 
scribed by a standard quantum optical master equation 



: {Hcsp - 



pHln) 



r[p], 



(2) 



where for our purposes we explicitly take into account the 
indistinguishability and the motional degrees of freedom 
of the atoms. The terms appearing in the effective non- 
hermitian Hamiltonian TJoff and in the recycling term T[p\ 
are discussed in this subsection. 

The effective Hamiltonian can be written as the sum of 
a hermitian and a non-hermitian part, i/eff = ^cff ~^-^cS ■ 
In 

^iff =^fll^-n (4„Cg„ -I- 4nCen) 
n 

mnm'n' 

the first term accounts for the motion of the atoms in the 
harmonic oscillator potential. Here c^n (cl.^) are the an- 
nihilation (creation) operators for an atom in the internal 
state /? e {e,^} and in the vibrational mode \n) of the 
harmonic oscillator. These operators obey the fcrmionic 



anti-commutation relations {cl^, cL^^i} = {cjSn, c^'n'} 

and {cBniCj3'n'} — Sfj^fj'Sn.n' ■ The second term in 
Eq. ([3]), a two particle operator, is the dipole-dipole in- 
teraction between the atoms, which is induced by the 
collective coupling to the radiation field and has to be 
taken into account when the distance between the atoms 
is comparable to or smaller than Aq. It is given by 
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where F is the single particle decay rate from |e) to \g) as 
defined above, ipm{x) are the position space harmonic os- 
cillator eigenf unctions and d^g denotes the dipole matrix 
element unit vector of the transition. The representation 
Eq. ([5]) for the function G{^) is valid in the hmit |^| <C 1 



which is fulfilled in the Lamb-Dicke regime. For broader 
traps with interatomic distances fco|a; — a;'| ~ 1 there are 
corrections to G(^) of order 1/^^ [3D]. 

The non-hermitian part of the effective Hamiltonian 
reads 
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= /dr!^iv(fc)i?;,^(fc)i?^,„(fe), (7) 



where 
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is the angular distribution of dipole radiation, k — 
{ki,k2,kz) denotes the unit vector of the photon 
emission direction, Rmn{k) = 11^=1.2,3 -^mjnjC^i) = 
Y\, (?7ij|e~' ^'^^^°'^'^°'i>\nj) the recoil matrix elements (see 
Appendix \m , tjj = fco yJh/{2Mvj) the Lamb-Dicke pa- 
rameter corresponding to the j'-axis, and aj [a,) the an- 
nihilation (creation) operators for vibrational excitations 
along the j-axis. The effective Hamiltonian Eq. (Is]) and 
([6]) is valid under the assumption that the system size 
is small compared to the distance which light travels on 
typical timescales of the system dynamics. This is always 
fulfilled for optical transitions and a system size on the 
order of an optical wavelength. 
The recycling term is given by 

mnm' n' 

The coefficients Rn'm,'m,n account for the momentum re- 
coil of the emitted photon during a quantum jump. 

For our initial state of interest with two atoms, both 
residing in the ground state of the three dimensional har- 
monic oscillator and one being internally excited, the 
density operator is given by po = CggCgo I vac) (vac| CeoCgo- 
In the recycling term for this initial state r[po], quantum 
jump terms that correspond to a decay of the excited 
atom while remaining in the ground state of the harmonic 
oscillator (i.e. terms with m — Q, m! = 0) are identically 

zero due to the fermionic properties cJ^qSo ~ ^go^go ~ ^^ 
reflecting the Pauli-blocking of this particular channel. 
The total initial decay rate Fos of the excited atom is 
given by the sum over all available decay channels 



Teff = E (30;5"^|r[Po]|ffO;5m) 
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where \gQ;gm) — c)^c)g^ |vac). Equation ^ is the gen- 
eralization of Eq. pi) to a 3D setup. In the Lamb-Dicke 



limit rjj <Si 1, to second order in rjj only decay to the three 
oscillator states with one motional quantum {\gG\gm) 
with \fn\ — 1) contributes to F^ft. Expanding the re- 
coil matrix elements in the small parameters rjj yields 
|i?io(fc)P = kJTfj+Oi'q'l) for the j-axis. The total initial 

decay rate in this hmit is Fcff — X]j=i 2 3 ^jV]^ + ^('7j )j 
where Uj are numerical coefficients of order 1 and de- 
pend on the relative orientation of the transition dipole 
matrix element to the j-axis. More details on these coef- 
ficients and on F^ff beyond the Lamb-Dicke limit can be 
found in appendix [Xj For the special case of an isotropic 
harmonic oscillator (i.e. i^i = 1^2 = ^^3) the particularly 
simple result Fcs = rj^T + 0{r]*) is found. 

Let us now return to the discussion of the collective ef- 
fects introduced by the dipole-dipole interaction (Eq. Q) 
and the cross-damping terms (Eq. (|6|), still concentrat- 
ing on the case of two atoms trapped in the Lamb-Dicke 
regime. The dipole-dipole interaction couples different 
vibrational levels of two-particle states where one atom 
is internally excited and the other is in its internal ground 
state. In evaluating the matrix elements in Eq. m, di- 
vergences at short relative distances appear, which reflect 
the fact that, in principle, one would have to solve the 
exact two-atom problem to capture the right behaviour. 
Effects of the short-range potential are important at dis- 
tances comparable to a few tens of Bohr radii |32] • How- 
ever, in order to estimate the order of magnitude for 
the leading terms it is sufficient to treat the atoms as 
independent particles being described by their internal 
state and center of mass wavefunction. The divergence 
of the dipole-dipole matrix elements can then be avoided 
in a standard way by introducing a spatial cutoff for the 
relative distance between the atoms, which yields a fi- 
nite result independent of the cutoff. By this procedure, 
the leading terms of the dipole-dipole matrix elements 
are found to be on the order of F/(100r7'^). Note that 
the classical scaling of the dipole-dipole interaction with 
the inverse distance cubed is refiected in the 77 depen- 
dence. Another process in this case leading to dissipa- 
tive redistribution between excited oscillator states are 
non-diagonal cross-damping terms contained in the non- 
hermitian Hamiltonian Eq. ([6]). The leading terms are 
of order r]^T. 

In the case of an experimental realization with alkali 
atoms where the linewidth F is large, the dipole-dipole 
interaction is the dominant source of imperfections and 
competes with the Pauli-blocking effect as it couples the 
motional state po to motional states that are not blocked. 
A compromise has to be found for 77: smaller values of 
77 lead to a larger blocking effect, but also to a larger 
dipole-dipole interaction. When estimating the rate at 
which the system leaves the initial Pauli-blocked state 
due to dipole-dipole interaction one must take care of 
divergences in the sum over all intermediate states. In 
situations where this rate is on the order of the original F 
even for large 77, the Pauli-blocking of spontaneous emis- 
sions would no longer be observable because the atoms 
will be transferred to a configuration of motional states 



that is unblocked on timescales comparable to the origi- 
nal spontaneous emission rate. 

In contrast, for an implementation with alkaline-earth 
atoms the small natural linewidth of the metastable ex- 
cited state makes the dipole-dipole interaction and the 
cross-damping processes negligible. Additionally, as will 
be argued in the next subsection, the impact of these 
effects on the spontaneous emission process which is in- 
duced by admixing an excited state with a large decay 
rate can be avoided by coupling far off-resonantly to this 
intermediate state. 



B. Alkaline earth- like atoms 

In this subsection we further elaborate on the experi- 
mental realization of Pauli-blocked spontaneous emission 
with alkaline earth-like atoms, commenting in more de- 
tail on the initial state preparation in a "magic" wave- 
length optical lattice, quenching of the metastable ex- 
cited state to induce spontaneous emission and the ap- 
plication of an external magnetic field to decouple the 
nuclear and electronic spin degrees of freedom. To be 
specific we present the experimental scheme choosing the 
isotope ^^^Yb, which due to its nuclear spin / = 1/2 ex- 
hibits a particularly simple level structure. The discus- 
sion, though, can be readily adapted to other fermionic 
alkaline earth(-like) species, such as ^^Sr. 

In Fig. [2] the relevant electronic states and experimen- 
tal steps are illustrated. The doubly forbidden "clock" 
transition (6s^)^So ^ (6s6p)'^Po is weakly dipole-allowed 
due to hyperfine mixing with higher lying P-states and 
has a natural linewidth 7/(277) ~ lOmHz I33l|34]. The 
nuclear spin decouples from the electronic degrees of free- 
dom for the ground state with total electronic angular 
momentum J = 0, and the two magnetic sublevels are 
given by \^So;mj — ±1/2). For the metastable excited 
state ^Pq the total angular momentum projection mp 
rather than the nuclear spin projection mj is a good 
quantum number due to the small hyperfine admixture 
of electronic angular momentum. Nevertheless, the mag- 
netic sublevels are still almost pure m/ eigenstates. The 
states with mj = ±1/2 are denoted as |t) and jj,), re- 
spectively. 

We choose the ground state \g f) = |^So;to/ = 1/2) 
and the metastable excited state |e f) = pPo! "^/ = 1/2) 
as the starting point from which the Pauli-blocking ef- 
fect can be observed (see Fig. [TJa)), since the tiny nat- 
ural linewidth 7 of the "clock" transition renders both 
direct radiative decay and dipole-dipole interaction be- 
tween these two states negligible. 

The requirement of equal trapping potentials for the 
states |e f) and \g f) can be met by operating the opti- 
cal lattice at the "magic" wavelength of 759 nm [331 ES] . 
At this wavelength, the AC polarizability and, thus, the 
light shift caused by the lattice laser is the same for the 
^Sq and the '^Pq state to first order in the laser intensity. 
In the regime of a deep optical lattice where tunneling 



between the sites is negligible, the spatial potential on 
each lattice site is well approximated by a 3D harmonic 
oscillator with trapping frequencies u. Vibrational fre- 
quencies of z^/(27r) = 90 kHz have already been realized 
in a ID optical lattice [3S]. 

Let us now turn to the discussion of state preparation 
and the required sequence of laser pulses. 

(I) Loading of the optical lattice with two atoms per 
site can be achieved by adiabatically ramping up the lat- 
tice potential on a cloud of ultracold atoms [371 ISH] ■ For 
a temperature lower than the motional level spacing, a 
band insulator with two atoms per lattice site in the mo- 
tional ground state |0) and with opposite nuclear spin 
projections |t) and \\.) will form. The resulting state on 
each lattice site is cj^-^oSio K^*^) (^^^ Fig. 2 I). 

(II) To prepare a suitable initial state in order to ob- 
serve Pauli-blocking, the atom in \g i) |0) is excited with 
a laser 7r-pulse on the "clock" transition with Rabi fre- 
quency f2. Recoil-free excitation is possible in the eas- 
ily accessible regime i^max ^ ^ ^ j where most con- 
veniently the laser direction is chosen along the axis of 
strongest confinement. In this limit, spontaneous emis- 
sion can be neglected and the vibrational sidebands are 
well resolved. Choosing (t+ polarized light flips the nu- 
clear spin due to hyperfine interaction and furthermore 
ensures that the other atom is not affected by the pulse. 
By tuning the laser on the carrier transition, the initial 
state of interest c^^qcLq |vac) is prepared. 

(III) Quenching the excited state offers the possibility 
to induce a tunable effective decay rate from |e f) to \g t) 
that exceeds the natural linewidth 7, and also to produce 
a decay that will not flip the nuclear spin state (which 
is not the case for intrinsic decay of the ^Pq manifold, 
as the transition is weakly allowed due to hyperfine cou- 
pling). A dressing laser couples the metastable excited 
state |e t) on the weakly magnetic dipole-allowed transi- 
tion to the intermediate state (6s6p)^Pi [39], which has 
a large linewidth rip/(27r) = 29 MHz. Decay from ^Pi 
to states other than the ^Sq ground state occurs with a 
negligible probability. 

One has to ensure that during the decay induced by 
this dressing the nuclear spin of the initially excited atom 
is not flipped due to hyperfine interaction, as the decay 
channel |e f) —> Iff i) is not Pauli-blocked. This can be 
achieved by dressing |e f) with a state in the ^Pi manifold 
that only contains the |t) component, i.e. a product state 
of the form \mj; mj — 1/2). The only decay channel for 
such a state is into \g t) as follows from the selection rule 
Amj = for electric dipole transitions. In the following 
we will discuss how such coupling to an |t) state can be 
accomplished by either (i) using a^ polarized dressing 
light to couple |e f) exclusively to the magnetic sublevel 
with maximal mp which is a mj = 1/2 eigenstate, or 
alternatively (ii) using ir or a~ polarized dressing light 
while decoupling nuclear and electronic spin in the ^Pi 
manifold with an external magnetic field in the Paschen- 
Back regime [H] . To this end we diagonalize the Hamil- 
tonian governing the ^Pi subspace including hyperfine- 
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FIG. 4. Zeeman diagram of the ^Pi manifold. The eigen- 
states approach product states of electronic and nuclear spin 
for a large magnetic field. The states suitable for dressing 
the metastable excited state |e t) without flipping the nuclear 
spin during the decay are |mj) with niF ~ 3/2, 1/2, —1/2. In 
the inset, the probability for decaying without nuclear spin 
flip is plotted as a function of the external magnetic field. 



and Zeeman interaction 

H = AI- J/h^ + gjUBJ ■ B/h - giUNl ■ B/h. (10) 

Here gj [10] and gi [Hj are the electron and nuclear g 
factors, A |4^ is the magnetic hyperfine constant and B 
an external magnetic field. The eigenstates, which are 
denoted \'mp)^ approach product states of electronic and 
nuclear spin for large values of B in the Paschen-Back 
regime ^bB 3> A as illustrated in the Zeeman diagram 
in Fig. |4] The sign ± labels the nuclear spin projection 
in this limit, i.e. \mp) — >■ |mj=mFTl/2;TO/=±l/2). 
One can couple to the three states \rn^) with m^ — 
3/2, 1/2, —1/2, which are suitable as intermediate states 
for the quenching process, with a"'", tt and a~ polar- 
ized laser light, respectively. Their expansion in terms 
of product states \mj]mi) is 

I™/?) = cl^p l"^J =mF - 1/2; mi = 1/2) 






- CX, „ \mj = nip + 1/2; m/ = -1/2) 



The i?-dependence of the coefficients is given by 
1 



A(i) - 



1 



X + ^mp 



;2 + ^mpx + 1 



where x — 2{gj^B + g//iAr)i3/(3|A|) is a dimcnsionless 
variable for the magnetic field strength. The probabil- 
ity for decaying from \m'^) to \g f) or \g |) is |c^j^P and 
jcjl^j^p, respectively (see inset Fig. W|. Note that (i) if 
one chooses to work with a cr"^ polarized dressing laser, 
a weak magnetic field defining the quantization axis is 



sufficient because |cj/2l^ is identically 1 for arbitrary B- 

values, and a^ polarized light couples |e f) exclusively to 
|3/2+) in the ^Pi manifold. On the other hand, (ii) for 
■K and a~ polarized light the external magnetic field has 



to be sufficiently strong so that |c 



±l/2l 



approaches unity 



and the Zeeman splitting becomes large enough so that 
coupling to \m~pi is negligible. For a magnetic field of 
0.05 T the probability to decay without nuclear spin flip 
is > 95 %. The arguments presented here for -'^^-'^Yb are 
also applicable to other alkaline earth-like species which 
may have a nuclear spin I > 1/2. Method (i) requires 
that states with maximal or minimal mi and mp values 
are used in the protocol whereas method (ii) is suitable 
for arbitrary m/ sublevels. The decoupling of nuclear 
and electronic spin in the limit of large B also holds for 
isotopes with / > 1/2 where additional quadrupole ef- 
fects have to be taken into account in the Hamiltonian 
Eq. (pl [29l[33]. 



After having discussed the required experimental steps 
let us turn to the resulting Pauli-blocked decay dynam- 
ics. First we show that this dynamics is described by a 
coupled set of effective rate equations for the populations 
in |e) and \g). Subsequently we provide a suitable set of 
experimental parameters. 

We are interested in a regime where the detuning of 
the dressing laser Adr is large compared to all the other 
characteristic frequencies which govern the system dy- 
namics: the linewidth Fip, the dipole-dipole interaction 
for two atoms in ^Sq and ^Pi, respectively, the dressing 
laser Rabi frequency il^^ and the trapping frequency ly of 
both the metastable excited and the ground state. Un- 
der these conditions the intermediate state is only virtu- 
ally populated and can be adiabatically eliminated. Fur- 
thermore, in the resulting effective decay dynamics the 
ground state atom in \g t) |0) acts as a spectator not 
taking part in the dynamics but blocking the dominant 
decay channel |e f) |0) -> \g t) |0) for the excited atom. 
The corresponding dynamics in the Lamb-Dicke limit and 
for the special case of an isotropic trap is described by 
the rate equations 

Pet = - (14. r iv' + vl) + 14. r) rPet + o{r,^) 



P. 






gi = l^m.l rPet- 



Here we have introduced the effective decay rate P ~ 
i7^j./(4A^j.)rip and the populations of the reduced sys- 
tem density matrix 

Pet = (vac \cetoCg^opc\^QC^^^Q\ vac) 
At= Y^ (vac|cg^m,Cgtop4.^o4tmKac) 

Pgi = Y^ (vac \cgim,Cg^opcl^ocl^J\ vac) . 

m 

The latter describe the probability of finding the initially 
excited atom in |e t)7 \g t) and \g |), respectively, while 



the other atom remains in |(7t)|0)- The two Lamb- 
Dicke parameters rj and r]di correspond to the momen- 
tum recoils of the induced spontaneous emission and the 
dressing laser photon, respectively. The main result for 
the initial state c^^gCg-KQ |vac) described by P^^ = 1 and 
Pg^ = Pgi = is a total decay rate from |e f) |0) given 
byT 



off 



= (\rt 



'(.V + Vdi) + Icmpl )r- In the regime of 
interest jcJ^^P — ^ 1 and jc^^P — > 0, the decay rate Toff 
is of order 77^ smaller than the effective induced decay 
rate F from |e f) to \g f) for a single atom. 

In this discussion so far we have not included elastic 
and inelastic collisions between the two atoms. Colli- 
sional shifts do not affect the experiment as long as they 
are small compared to the motional level spacing. State 
changing collisions can lead to loss of the atoms from 
the lattice. While, in principle, the coUisional loss from 
two atoms in ^Pq could be strong, this situation is never 
encountered in our protocol. Here, we only require the 
^Sq-'^Po inelastic collisions to be small compared to the 
quenching rate in the experiment. 

A suitable choice for the detuning Adr is dictated by 
the requirement that this has to be the largest frequency 
scale in the dynamics and simultaneously a sufficiently 
large effective decay rate F has to be induced for a given 
dressing laser Rabi frequency ftdr- For ^^^Yb with a typi- 
cal trapping frequency i'/{2tt) = 90 kHz, the Lamb-Dicke 
parameters rj = 0.28 and rydr = 0.09 correspond to tran- 
sitions with wavelengths A = 399 nm and Adr = 1285 nm, 
respectively. The dipole-dipole interaction between two 
atoms in ^Sg and the virtually populated ^Pq, respec- 
tively, can be estimated to be on the order of Fip so 
that Adr 3> Fip is required. For Adr = lOFip, with a 
Rabi frequency ildr = 4 MHz, an effective decay rate of 
F = 220 Hz can be reached. Note that as an alterna- 
tive to the direct coupling ^Pq^-j-^Po, quenching of the 
metastable excited state can be done via a two-photon 
Raman process involving the dipole-allowed transition 
■^PqO'^Si and the intercombination transition '^SiO^Pi 
which can lead to a larger two-photon Rabi frequency 

Odr ESI. 



C. Alkali atoms 

Let us now turn to alkali atoms with their fastly de- 
caying low lying transitions (see Fig. 3]). Here, both the 
laser excitation and the subsequent decay dynamics are 
fast compared to the motion of the atoms in the trap. In 
contrast to the alkaline earth case, where it is important 
to have equal trapping potentials for ground and excited 
state atoms, the notion of motional eigenstates becomes 
irrelevant for the spectrally very broad excited states in 
alkali atoms, therefore the optical potential for internally 
excited atoms does not play a crucial role in this scenario. 

First we will analyze the experimental requirements 
and sequences to observe Pauli-blockcd spontaneous 
emission using alkali atoms. Subsequently, we will dis- 
cuss imperfections associated with the large linewidth of 



the considered transition and estimate their impact. 

The relevant level scheme and the sequence of laser 
pulses is sketched in Fig. [s] for '^''K as a represen- 
tative of the alkali family, though the scheme pre- 
sented here is not specific to a particular species and 
can readily be adapted to other fermionic alkali iso- 
topes. We propose to use the magnetic sublevels 
with maximal projection of total angular momentum 
niF, i.e. \g) = pSi/2; i^ = 9/2; wp- = 9/2) and |e) = 
pP3/2;-P = 11/2; m_F = 11/2), as ground and excited 
states as depicted in Fig. W[h) ■ For atoms in the inter- 
nal state |e), the decay to\g) is the only dipole-allowed 
transition, ensuring that these states form a closed ef- 
fective two-level system regarding spontaneous emission. 
To prepare an initial state suitable to observe the Pauli- 
blocking effect, the following experimental sequence can 
be used. 

(I) The lattice is adiabatically ramped up on an ul- 
tracold cloud of atoms which are internally in a 50% 
mixture of the state \g) and the magnetic ground state 
sublevel |i) = pSi/2; F = 7/2; mp = 7/2), which leads to 
a band insulator with two atoms per site in the internal 
states \g) and \i), respectively, provided the temperature 
is lower then the motional level spacing [37l [38] . For suf- 
ficiently deep lattices, when tunneling is negligible on the 
experimental timescale, the state on each lattice site is 
l''/'i) — 4o'^gO |vac). (II) A TT-pulse is applied on the tran- 
sition \g) — ^ |e) with a Rabi frequency fli to internally 
excite one of the atoms. On the one hand, we require 
r2i 3> F so that spontaneous emission during the excita- 
tion can be neglected to lowest order. On the other hand, 
fli has to be small compared to the hyperfine splitting 
between the states \g) and |i), so that the other atom 
is not affected by the pulse. As the ground state hy- 
perfine splitting typically is on the order of GHz, both 
requirements can be fulfilled simultaneously. For alkali 
atoms r^i ^ F also implies fii 3> f, thus the trans- 
fer takes place in the strong-excitation regime [21]. The 
initial motional wavepacket is essentially unchanged dur- 
ing the excitation, apart from a momentum kick hk]^ by 
the absorbed photon, where k^ is the laser wave vec- 
tor. The state prepared by the laser pulse is therefore 
\ipu) = X^Ti^^'^Io'-en KS'C) where the expansion coefR- 
cients r^ are related to the recoil momentum associated 
with the absorption of a laser photon: r„ ~ (nje'^'^'^IO). 
For tight traps, in the Lamb-Dickc limit, the probability 
to find the excited atom in its motional state |0) after 
the laser pulse is |rop ~ 1 — ^ . kj^^rj'j in leading order 
in r]j. (HI) A second 7r-pulse transfers the other atom 
from \i) to \g) with a radio frequency pulse or by a two- 
photon Raman transition with effective Rabi frequency 
^2 3> F. Spontaneous emission from the excited state is 
neglibigle during the whole pulse sequence when the total 
time needed to perform both pulses is short compared to 
F^^. As the energy gap between the initial and final state 
for the second pulse is not an optical frequency but the 
hyperfine splitting of the ground state, momentum recoil 



can be neglected in this case. The state prepared after 



the second pulse is IV'iii) = J2 



n n'-gO'-en 



TnC. 



I vac). This is 



the initial state for the Pauli-blocked decay dynamics. 

The time evolution of the initial system density matrix 
Pa ~ IV'iii) ("^iiiI is determined by the master equation 
Eq. (pi), for both shallow or tight confinement in the lat- 
tice. The initial effective decay rate can be calculated in 
analogy to Eq. (lo]), with IV'iii) as the initial state. For 
the special case of an isotropic trap in the Lamb-Dicke 
limit it is found to be TeS = 2t]^T + 0(77"*), where the 
factor 2 reflects the fact that there is a momentum kick 
involved both during the initial state preparation of IV'iii) 
and during the spontaneous emission of a photon. 

The dispersion of the excited atom's motional wave 
packet due to the different optical potential for the 
excited state can be described by an additional term 



Z-,n,n' "'-^n'.n'-^en' ' 

The rates T^, „ 



' " in the effective Hamiltonian Eq. (|3| . 
are typically on the order of v and can 



be neglected compared to the effective decay rate. 

Let us return to collective effects in the master equa- 
tion. The dipole-dipole interaction induces transitions 
away from the initial state towards motional states which 
are not Pauli-blocked and decay essentially with the rate 
r. Therefore, it competes with the effect one wishes 
to observe. A compromise has to be found regarding 
the confinement strength: tighter confinement and thus 
smaller values of 77 are preferable on the one hand in order 
to suppress recoil induced change of the initial motional 
wavepacket, but lead to more pronounced dipole-dipole 
interaction on the other hand. For the alkali scheme 
where the atoms are prepared in internal states which en- 
counter strong dipole-dipole interaction, the rate of tran- 
sitions of the system away from Pauli-blocked states is 
significant. Estimating these rates is difficult because of 
diverging matrix elements which emerge if the interaction 
potential at short interatomic distances is not properly 
treated. For cutoffs of the order of several hundred Bohr 
radii, which is a typical value for Potassium, simple esti- 
mates show that the timescale for transitions away from 
Pauli-blocked states is on the order of the original F even 
for 77 of order one, with shorter cutoffs leading to larger 
values. To determine whether this is realistic in experi- 
ments would require a much more detailed treatment of 
the interatomic potential at short scales. Avoiding real 
population of states which encounter a large dipole-dipole 
interaction is an advantage of the alkaline-earth scheme 
compared to the alkali scheme. 



IV. SHAPE OF THE EMITTED PHOTON 
WAVEPACKET 

In the previous sections we have discussed the atomic 
dynamics of the Pauli-blocked decay. In addition, this 
effect also becomes manifest in the intensity distribution 
of the emitted photon wavepacket. Compared to the de- 
cay in absence of a blocking atom, due to the prolonged 
lifetime of the excited state, we will see a broader spatial 




FIG. 5. (a) Energy diagram of the initial state: two fermions 
with blocking atom in a superposition of the motional states 
jo) and jl). (b) Spatial intensity distribution It{x) of the 
emitted photon wavepacket at time t for the limiting cases 
where the blocking atom is in |0) or |1), respectively. In both 
cases, the outgoing wavepacket propagates with the speed of 
light and falls off exponentially with a width c/((l — ari^)r) 
(c/(Q?7^r)) for the initial state |3l;eO) (|c/0;eO)). (c) The 
blocking atom is prepared in a superposition of motional 
states, we choose fig = ^ ~ if /"^ B.nA /ii — 1^/1 — fJ-o- As 
this is the motional state which the excited atom would reach 
by sending out a photon along the negative i-axis, emission 
in this direction is suppressed at the beginning of the atomic 
decay (i = 0). Hence, the photon intensity is zero at a: = —ct 
and has a maximum at 2; = +ct. During the decay of the 
excited atom, the blocking atom oscillates in the trap with 
a frequency u, which imprints the spatial period 2Trc/u onto 
the photon wavepacket. 



extension of the wavepacket. Furthermore, to some ex- 
tent a temporal shaping of the emitted photon is possible 
by preparing the blocking atom initially in a superposi- 
tion of motional states. For simplicity, we regard the 
case where the blocking atom is in a superposition of the 
lowest two motional states (see Fig. Isja)). This situation 
captures the essential effect as the contribution of higher 
motional states in the superposition leads to modifica- 
tions of order 77^ or less. The photon shaping effect can 
be most explicitly observed in the regime T <^ v, which 
can be realized in the setup with alkaline earth atoms as 
discussed above. In this scheme, the appropriate initial 
state could be prepared after step II, by coherently trans- 
ferring part of the population of the blocking atom from 
the motional state |0) to |1). 

We study this model system in a Weisskopf-Wigner 
ansatz, where we have one excitation either in the 
atomic system or in the radiation field as an one-photon 
wavepacket. Compared to the previous treatment where 
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we have traced over the radiation field, here the infor- 
mation about the photon is contained. We consider mo- 
tional excitations of the atoms only in one dimension and 
additionally, for the excited atom, only include the low- 
est motional band |0) as shown in Fig. Isja). Therefore, 
we do not include the effects of the dipole-dipole interac- 
tion and cross-damping. As discussed above, in the limit 
r <C i^, which corresponds to the alkaline-earth case, the 
effect of these interactions can be suppressed by choosing 
a sufficiently large detuning in the quenching process. 
Within this ansatz, the state can be written as 

IV-W) = "oe W IffO; eO) -I- aie(t) I5I; eO) 

+ ^amn,fc(i) |ff'7i;5n;lfc) , (11) 

n>Tn,/c 

where the first two terms correspond to the situation be- 
fore the decay and the last term corresponds to both 
atoms in their internal groundstate and an outgoing one- 
photon wavepacket. The initial condition is aoe(O) ~ 
Ho, aie(O) = /xi, afnn,k{0) — for all m,n,k and 
ImoP + |y"iP = 1- We solve the Schrodinger equation 
of the Hamiltonian 



H 



y jvjjQi 






cre{e,g},n 
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fihJkb\hk 



E 



8kd^ 



^kbk > Rn 



n^m 



lik-xH 



which consists of internal and center of mass degrees of 
freedom of the atoms, the free radiation field and the 
coupling between atoms and field in dipole and rotating 
wave approximation. Here, in addition to the operators 
and parameters defined before, bk (bk) is the annihila- 
tion (creation) operator for a photon with frequency Wfc 
and wave vector (polarization) k (e^), x is a unit vector 
in direction of the ID harmonic oscillator potential and 
£k — iy^hojk/(2eoL^) with quantization volume L'^. The 
solution for the coefficients aoe(t), aie{t) and amn,k{t) is 
obtained by the standard resolvent method [26 and is 
given in Appendix [B| 

From the spectral distribution of the photon amn,k{t) 
we calculate the first order correlation function I{r,t) — 
{tp{t)\E^-'> {r)E'^+\r)\ip{t)) , which is proportional to the 
electric field intensity. Here, E^-^^r) = ^^^k^k^^bk^^^'^ 
and its hermitian conjugate E'^'\r) is the positive and 
negative part of the electric field operator, respectively. 
The probability of detecting a photon between the times 
t and t + dt in a volume element r^dr dfJ around r is 
proportional to /(r,f)r^dr dfJdi [15]. 

Let us consider the atomic system placed at the origin 
of our coordinate system, with d^g parallel to the £-axis 
and X parallel to x. We analyze I{r,t) in the far zone 
limit kQr 3> 1 and for times t 3> T~g , when the initially 
excited atom certainly has decayed to the ground state. 
In the following, we discuss the properties of /(r, t) along 
the i-axis where we denote I{r,t)r'^\r={x.o,o) = It[x). 
The complete form of I{r,t) can be found in Appendix 

E 



For the two limiting cases of the initial state, where all 
of the population of the blocking atom is either in the 
motional state |0) or |1) (see Fig.lsjb)), the decay of the 
excited atom will be maximally blocked or essentially not 
blocked, respectively. The corresponding decay rates to 
order 0{tj^) are r|,^' = arj^T and F^J^^ == {l~ar]^)r. Con- 
sequently, the outgoing exponential photon wave packet, 
which has a spatial width c/F^jj' and propagates with 
the speed of light c, is much broader for the initial state 
\gO; eO) compared to the case of an initial state \gl; eO). 

Let us now consider the case where the blocking atom 
initially is in a superposition of motional states. As such 
a state is not an eigenstate of the harmonic trapping 
potential, the blocking atom oscillates in the trap while 
the decay of the excited atom takes place. Due to this 
osciallation, the preferred direction of the emitted pho- 
ton varies in time. The characteristic frequency for the 
dynamics in the trap is v, which is imprinted onto the 
photon wavepacket as a spatial oscillation with a period 
2ttc/i' (see Fig.lsljc)). 

We finally point out that other techniques to induce 
dynamics for the blocking atom could allow for more so- 
phisticated shaping of the emitted photon. For example, 
Rabi oscillations could be driven for the blocking atom 
between \gO) and another, decoupled internal state, or 
the trapping potential could be modulated during the 
Pauli-blocked decay. 



V. SUMMARY AND OUTLOOK 

In this work we have studied how the Pauli exclu- 
sion principle can give rise to a suppression of sponta- 
neous emission from electronically excited states for cold 
fermionic atoms stored in optical lattices. Here the pres- 
ence of both a ground state and an excited state atom 
at the same lattice site can block the dominant decay 
channel for the excited atom, thereby significantly de- 
creasing the spontaneous emission rate. Complementary 
to the atomic dynamics we have studied how the decay 
dynamics in the presence of a blocking atom manifests 
itself in the characteristics of the emitted photon. We 
have suggested and analyzed experimental realizations 
with alkaline earth atoms and also with alkali atoms. 

On the one hand, from a conceptual point of view, 
observing these effects for fermionic atoms in optical lat- 
tices would be the first experimental demonstration of 
Pauli-blocked spontaneous emission. On the other hand, 
from a more practical perspective, the use of this block- 
ing effect in a controllable way can constitute an addi- 
tional, valuable tool in the context of reservoir engineer- 
ing in cold atom systems. Here, the idea of engineering a 
controlled coupling to an environment for the dissipative 
preparation of entangled states and many-body quantum 
phases has been explored both theoretically [T7H22] and 
in experiments with atomic ensembles |23j and trapped 
ions |24j . In particular, for cold fermionic atoms in opti- 
cal lattices Diehl et al. PHI have proposed and analyzed 
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a scenario, where quasi-local single-fermion dissipative 
processes can be tailored such that they lead to "cool- 
ing" into a BCS-type state of d-wave symmetry. The 
corresponding set of fermionic quantum jump operators 
for the dissipative dynamics are suggested to be imple- 
mented via a stroboscopic sequence of coherent and dis- 
sipative steps in a system of two-component fermionic 
alkaline-earth atoms. Here, controlled induced sponta- 
neous decay, which is suppressed (enabled) in the pres- 
ence (absence) of a second fermionic atom, as studied in 
the present work, constitutes the dissipative ingredient, 
which also warrants the required Fermi statistics of the 
quantum jump operators. 

In addition, the possibility to control to some extent 
the spatial and temporal emission characteristics of the 
outgoing photon under induced spontaneous decay in the 
presence of a suitably prepared blocking atom, might con- 
stitute an interesting tool in the context of the develop- 
ment of (directed) single-photon sources [ii] . 



where 
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x^ 



is the generalized Laguerre polynomial. For a tight trap- 
ping in the Lamb-Dicke limit, the recoil matrix elements 
can be expanded in the small parameter rjj <§; 1. To 
order rj^, the result is 



R 



\ — [m + \)k^rl^ for mj = Uj 

'ikjTijy'rn for \mj — nj\ — 1 

-^k'^r]'^ y/m,{m — 1) for |?7i,- — n,| = 2, 



2"'j'ij V ""V"- 1) foi' I™--' 

where we have defined m = max(?77,j , Uj 
matrix elements are those with m,- 






. The dominant 
whereas matrix 

elements with |mj — nj\ — 1 and \mj — nj\ = 2 are of 

order rjj and ry? respectively. 

The master equation coefficients Rn'm'mn defined in 
Eq. ([7| are real, because of parities they are zero if any of 
the three components oi n' + m' +Tn+n is odd, and they 
are invariant under exchange [n'm') <->■ irnn) and under 
exchange n' -(->■ ?7i'- or rrij -H- Uj for any j' = 1,2,3. In 
Lamb-Dicke expansion to order 77^, the only contributing 

m = n and m' = n', (ii) 
jj for an i G {1,2,3}, (iii) 
2(5i J- for an i E {1,2,3} (or vice 
In these cases, we find 



coefficients are those with (i 
IttIi — nA — \m': — n'A = 5. 



m = n and 



lm;.-n;.| = 



versa 



J2j{^'j +™j + l)"j'?i 



Rn'm'mn ~ { tttTjf Vmhn 

-\(^i'ni\/'m'{m' - 1) 



(i) 

(ii) (Al) 

(iii), 



where m! = max(?7i^,n^), m — Taax{mi,ni). The numer- 
ical coefficients aj = ^(2— MegjP) with | < a^ < | and 
^ aj — 1 depend on the projection of the transitions 
dipole matrix element onto the j-axis. 

The total initial decay rate Tes defined in Eq. ^ can 
be written as 



Appendix A: Recoil matrix elements and initial 
Pauli-blocked decay rate 



In Eqs. (|6| and (Is]), the recoil matrix elements 



Ecff — r (1 — -Roooo) , 

which in the general case can be evaluated numerically. 
For the case of tight trapping in the Lamb-Dicke limit, 



by using Eq. ( |AT| ), it is given by Feff = I]j=i,2,3 "j'7|r - 



o(r?l). 



-t^mjUj yf^j 



-^rl^/\ 



/ min(m„n,)! . |„^._„^.| 
max(mj, rij)'. 



Appendix B: Weisskopf-Wigner approach 



For times t ^ T 



off 



Wigner ansatz Eq. (Ill in 

aoi,fe(^), aon,fc(i), ain,fcT^ with n > 1 is 



the solution of the Weisskopf- 
terms of the coefficients 



12 



aoi.kit) — ^k^eg ■ e^e 



<10n,k{t) — ^kdeg ' ^k^ 



-i{u>u+v)t 



lioRm{k-x) 



^^lRmik■x) 






ai„,.(t) = 5,!de, • e,e-'(-'=+("+i)'^)* '^^^^-^'^ ' * 



/i (K + {n + l)u - wo) + ir^ffVs) ' 



r 



.(0) 



„4^ ^T-(i) 



where T^g' = ayy^r + ©(ry") (r^^^ = (1 - a772)r + 0{r^^)) (I5I; eO)). All other coefficients a,„„,fc(t) with n > m 
is the effective decay rate for the initial state jgO; eO) are zero within this model. The first order correlation 

function I{r,t) can be written as [45] 



r(r,t) = |*i(r,OP + E[l*i°Hr,t)P + |*i^nr,t)|^ 



n>l 

„ (R^u)^ sin^ 9 



(47reo)^c^r2 



-2 (^oi?Si(")Mii?oo(u) exp (iz.(i ~ r/c)) exp (-(r'"^^ + rW)(i _ r/c)/2)) } 



{iTreo)'^c'^r 



I 

where 9 is the angle between r and d^g , u = r ■ x and 
is the Heaviside step function. 
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